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The H+ and H 2 D+ ions are important probes of the physical and chemical condi- 
tions in regions of the interstellar medium where new stars are forming. This paper 
reviews how observations of these species and of heavier ions such as HCO + and 
HaO + can be used to derive chemical and kinematic properties of nearby pre-stcllar 
cores, and the cosmic-ray ionisation rate toward more distant regions of high-mass 
star formation. Future prospects in the field are outlined at the end. 
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1. Introduction 

The formation of stars and planets is a crucial step in the life cycle of matter in 
galaxies, but remains poorly understood. While optical and near-infrared observa- 
tions are useful to locate newborn stars and measure the distribution of ionised gas, 
the cool and dense gas where the actual star formation takes place is best probed 
by long-wavelength (infrared through radio) observations. Molecular spectroscopy 
is an essential part of this effort, thanks to the richness of molecular spectra with 
many rotational and rovibrational lines which probe a large range of temperatures, 
densities, and optical depths. 

This paper reviews the role of line observations of and its isotopologues in 
the determination of basic parameters of star-forming regions. After an outline of 
modern methods to obtain descriptions of the temperature and density structure 
of star- forming regions (§ 2), § 3 shows how observations of (and other ions) 
may be used to measure the ionisation rate of molecular gas. Section 4 discusses 
determinations of the chemistry and kinematics of pre-stellar cores from H2D + and 
D2H" 1 " observations, and § 5 concludes with an outline of future opportunities. 



2. Measuring molecular abundances 

The interpretation of molecular spectra from star-forming regions requires knowl- 
edge of their physical structure. The envelopes of young stars have temperature 
and density distributions which are strongly centrally peaked. For example, in the 
envelope of AFGL 2136, one of the first sources where Hjj" was discovered, the 
temperature drops from 450 K at 6xl0 15 cm to 22 K at 2xl0 18 cm from the cen- 
tral star, while the density drops from ~10 7 to ~10 4 cm~ 3 over the same distance. 
Such gradients lead to significant chemical stratification. In the cool outer envelope, 
the gas-phase chemistry is dominated by ion-molecule reactions which do not have 
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activation barriers. Closer to the central star, neutral-neutral reactions become in- 
creasingly important. In addition, at low temperatures, atoms and molecules freeze 
out on grain surfaces where they may react with each other. The dominant type of 
reaction is hydrogenation at low densities and oxygenation at higher densities. The 
resulting ice mantles evaporate off the dust grains when the temperature exceeds 
-100 K. 

The temperature and density gradients in protostellar envelopes also manifest 
themselves as gradients in the molecular excitation. While in the case of infrared 
observations of rovibrational bands, the excitation temperature can be measured 
directly, the determination of molecular column densities from radio observations 
often relics on a single or just a few transitions. To be able to compare column 
densities from different molecules which may have been derived from lines with very 
different upper level energies and/or radiative decay rates (hence critical densities), 
a comprehensive excitation model of the source is essential. 

Modern schemes to construct models of star-forming regions are discussed by 
Doty et al. (2004). Total mass and temperature structure of the envelope is obtained 
from dust continuum observations, while the density structure is inferred from 
molecular lines. In constraining model parameters, it is important to be aware 
that different lines are observed with different beam sizes and thus do not sample 
quite the same gas. For further details please consult Van der Tak (2005). 

3. Probing the cosmic-ray ionisation rate 

One key parameter of star-forming regions is the ionisation fraction of the molecular 
gas. First, the electron fraction determines the influence of magnetic fields on the 
dynamics, in particular the efficiency of magnetic support against gravitational 
collapse. Furthermore, the ionisation fraction sets the time scale of ion-molecule 
chemistry, which is the most important type of chemical reactions in the gas phase 
at temperatures of — 10K and densities of — 10 4 cm~ 3 . Under these conditions, the 
bulk of the ionisation is due to 'cosmic rays' (—100 McV protons, produced in 
supcrnovae), at a rate (cr, while ionisation by stellar— 1-10 keV X-rays becomes 
important on small scales (Stauber et al. 2005). To estimate £cr, Van der Tak & van 
Dishoeck (2000) have used models of the temperature and density structure of the 
envelopes of seven young high-mass stars to model observations of Hg absorption 
and HCO + emission lines. They find £cr ~ 3xl0~ 17 s -1 , in good agreement with 
measurements of the local flux of low-energy cosmic rays with the Voyager and 
Pioneer spacecraft (Webber 1998). The models reproduce —half of the observed 
V(Hg ), suggesting that the other half arises in intervening clouds. 

Recent observations suggest strong variations of £cr between lines of sight, 
which appear related to column density. On the one hand, observations of (sub-)- 
millimctrc emission lines of DCO + and other ions towards the low-mass pre-stcllar 
core LDN 1544 indicate a decrease of £cr by a factor of —10 (Caselli et al. 2002), 
although depletion of CO and other neutrals onto grains may influence this result. 
On the other hand, observations of Hg absorption lines toward the £ Per diffuse 
cloud (McCall et al. 2003) indicate an enhancement of £cr by a factor of —10 
(Le Petit et al. 2004). These estimates are quite reliable now there is a consensus 
between experiment and theory on the rate coefficient of dissociative recombination 
of H^ at low temperatures (see Greene, this conference). Even larger cosmic-ray 
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ionisation rates of several 10 5 s 1 are found towards the Galactic Centre, where 
Oka et al. (2005) obtain 7V(H+)-4xl0 15 cm 2 from multi-line H^ observations 
(see also Geballe, these proceedings). However, this TN^H^) estimate may change 
slightly if infrared pumping is taken into account. The high stellar density in the 
Galactic Centre implies an infrared radiation field which is ^1000 times higher 
than in local clouds. Furthermore, the assumption of steady-state excitation of HjJ~ 
may break down at the high Hjj~ formation and destruction rates in the Galactic 
Centre clouds. Calculations of N(lt^ ) and (qr toward Sgr A which take these 
effects into account are needed (see also Le Petit, these proceedings). Nevertheless, 
an enhanced ionisation rate in the inner 250 pc of our Galaxy is consistent with its 
strong synchrotron radio emission which signifies an enhanced electron component 
of the cosmic rays. Furthermore, the Galactic Centre is a strong source of X-rays 
which may provide additional ionisation. 

Even without the Galactic Centre region, evidence is mounting that £cr is 
higher in regions with lower column density. Cosmic rays penetrate to a typical 
surface density of ^50gcm~ 2 (Umebayashi & Nakano 1981), which corresponds to 
an iV(H 2 ) of ^10 25 cm~ 2 . Even the Sgr B2 clouds, which arc the most massive 
clouds in the Galaxy, only have A^H^) ~ l-3xl0 24 cm~ 2 , so that absorption of cos- 
mic rays probably only plays a minor role except very locally. Instead, the apparent 
scaling of the ionisation rate with density may be due to cosmic-ray scattering off 
plasma waves. This process is more efficient in denser clouds where magnetic fields 
are stronger (Padoan & Scalo 2005), which would explain the observed trend. 

Liszt (2003) has suggested that the models for star-forming regions fail to in- 
clude recombination on dust grains. This process is certainly important for diffuse 
clouds, where grains are small and strongly positively charged by UV photons (Van 
Dishoeck & Black 1986). In dense star- forming regions, coagulation increases the 
average size of the grains and thus reduces their abundance by number, given the 
observed constant mass ratio of dust and gas (Ossenkopf & Henning 1994). How- 
ever, grain charging in these regions is mainly by collisions, and the negatively 
charged dust grains are very good at neutralizing positive gas-phase ions (Tielens 
2005). The processes governing the size and charge of the dust in dense clouds are 
thus quite different from those in diffuse clouds, and detailed models are needed to 
assess their importance for the overall ionisation balance. 

One critical parameter of chemical models of star-forming regions is the ratio of 
the He and H ionisation rates. Only the He ions are able to break strong chemical 
bonds such as the triple bond of the CO molecule, so that the atomic carbon 
abundance for instance depends sensitively on Cae/Ca- This ratio is insensitive to 
the slope of the cosmic-ray energy spectrum since most damage is done by secondary 
electrons, whose spectrum depends on the electron fraction and the H/H2 ratio of 
the gas. Observational constraints on these quantities are urgently needed. 

4. Chemistry and kinematics of pre-stellar cores 

Lacking an allowed rotational spectrum, the H3 ion needs a mid-infrared 'back- 
ground lamp' to be visible in absorption. Using H 2 D + as an alternative route to 
Ccr is problematic because the strong temperature dependence of the H 2 D + /H3 ra- 
tio means that the two species do not trace the same gas. Other ions such as HCO + 
and H 3 + are more suitable as substitutes of H3 , particularly if observations of CO 
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and H 2 arc available (Van der Tak et al. 2006). Instead, the astro-physical use of 
H 2 D + is that of a 'chemical filter' to image the distribution of cold ( 1 10 K) dense 
(> 10 5 cm~ 3 ) gas. Under these conditions, the formation of H 2 D + from HD + H3" 
is strongly enhanced because the back reaction has a barrier; and its destruction is 
reduced because the main destroyers of H 2 D+, CO, O and N 2 , freeze out onto the 
surfaces of dust grains. See Roberts & Millar (2000) and Roberts (these proceed- 
ings) for a detailed description of these processes. 

The detection of strong H 2 D+ 372 GHz emission toward the pre-stellar core 
LDN 1544 indicates the existence of a region of several 10 16 cm radius where all 
elements heavier than He are frozen out onto dust grains. Under these extreme con- 
ditions, H^ and its isotopologues become the dominant carriers of positive charge, 
a role normally played by metals (S + , Fe + ) and CNO molecules (HCO + , R30 + , 
N 2 H + ). Additional evidence for CNO-free zones comes from the detection of D 2 H + 
(sub-)millimetre emission toward the pre-stellar core LDN 1689N (Vastel et al. 
2004) and other cores (Vastel, these proceedings). However, observations of deuter- 
ated ammonia toward the same cores indicate that the depletion may not be as 
complete as thought before (Lis et al. 2006), or limited to a smaller region. In- 
deed, the H 2 D+ abundance of ~lxl0~ 9 in LDN 1544 by Caselli et al. (2003) has 
probably been overestimated by up to a factor of 10 through the collisional rate 
coefficient (Van der Tak et al. 2005) and the kinetic temperature (Crapsi et al., in 
prep.). Nevertheless, the H 2 D + abundance is clearly higher in pre-stellar cores than 
around low-mass protostars (Stark et al. 1999, 2004). 

Thus, although better data are needed to determine the sizes of the CNO-free 
zones, the very existence of such zones at the centres of pre-stellar cores implies 
that the relative abundances of H3 , H 2 D + and D 2 H+ are sensitive to the grain 
size (Flower et al. 2004) and the o/p H 2 ratio (Flower et al. 2006). Normally, these 
parameters are very hard to constrain by observations because of other competing 
effects. The centres of pre-stellar cores thus form a unique laboratory to study basic 
physical parameters of interstellar gas. 

In most radio observations, the lines are spectrally resolved, and the line profiles 
can be used to extract additional information about the source. A prime example is 
the double-peaked line profile of H 2 D+ toward the pre-stellar core LDN 1544 (Van 
der Tak et al. 2005). The central absorption dip implies a gradient in the excitation 
of H 2 D + , which means an increasing density and a constant gas temperature. Pre- 
vious models with a constant density at the centre and a centrally decreasing gas 
temperature are ruled out by the new H 2 D + observations. Furthermore, the nar- 
rowness of the absorption from the very light H 2 D + molecule implies that radial 
motions are subthermal in the outer parts of the core but increase inward. Maps 
of H 2 D + and D 2 H + emission will provide further constraints on the structure and 
kinematics of pre-stellar cores (Vastel, these proceedings) . 

5. Future directions 

The next 5-10 years will bring several new facilities for (sub-)millimctrc observa- 
tions of H 2 D + and D 2 H + . Specifically, the APEX telescope can observe the H 2 D + 
372 GHz line about ten times faster than the CSO (Belloche et al., in prep.). Its 
recently commissioned CONDOR receiver for Terahertz frequencies enables us to 
observe the 1370 GHz line (Emprechtinger & Wiedner, this conference). The D 2 H+ 
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molecule may be observed through its 692 GHz line with the CHAMP heterodyne 
array receiver, while its 1477 GHz line can be reached with CONDOR. Ultimately 
the ALMA interferometer will provide sensitive high-resolution images of the H 2 D + 
372 GHz and D 2 H+ 692 GHz lines. 

At mid-infrared wavelengths, higher spectral resolution observations (such as 
offered by CRIRES on the VLT) will be crucial to obtain velocity-resolved pro- 
files of the H^ absorption. Velocity information is crucial to distinguish intrinsic 
absorption from intervening clouds, and will improve the accuracy of estimates of 
the cosmic-ray ionisation rate. Observations of H 2 D + , D 2 H + and mid-infrared 
absorption lines would allow much more direct probes of the H 2 D + /H3 ratio. First 
attempts were unsuccessful but did result in the serendipitous infrared detection of 
formaldehyde (Roueff et al. 2006). 

Theoretical calculation of the rate coefficient of de-excitation of H 2 D + in colli- 
sions with H 2 would greatly improve the accuracy of estimates of the H 2 D + abun- 
dances in pre-stellar cores and star-forming regions. Such calculations are challeng- 
ing because of the reactive nature of the interaction, but the He-H 2 D+ system may 
be a useful starting point. 

In summary, H;j~ is great, but using it as an astrophysical tool requires other 
molecules too. The future will bring excellent observing facilities to carry out such 
observations. 

The author thanks Takeshi Oka for the invitation to this meeting, and John Black, Carsten 
Dominik, Eric Herbst, Harvey Liszt and Malcolm Walmsley for fruitful discussions. 
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